Coulomb blockade conductance peak fluctuations in quantum 
dots and the independent particle model 



Raul 0. Vallejos 1 , Caio H. Lewenkopf 1 , and Eduardo R. Mucciolo 2 

1 Instituto de Fisica, Universidade do Estado do Rio de Janeiro, 
R. Sao Francisco Xavier 524, 20559-900 Rio de Janeiro, Brazil 
2 Departamento de Fisica, Pontificia Universidade Catolica do Rio de Janeiro, 
Caixa Postal 38071, 22452-970 Rio de Janeiro, Brazil 
(February 1, 2008) 



Abstract 

We study the combined effect of finite temperature, underlying classical dy- 
namics, and deformations on the statistical properties of Coulomb blockade 
conductance peaks in quantum dots. These effects are considered in the con- 
text of the single-particle plus constant-interaction theory of the Coulomb 
blockade. We present numerical studies of two chaotic models, representative 
of different mean-field potentials: a parametric random Hamiltonian and the 
smooth stadium. In addition, we study conductance fluctuations for different 
integrable confining potentials. For temperatures smaller than the mean level 
spacing, our results indicate that the peak height distribution is nearly al- 
ways in good agreement with the available experimental data, irrespective of 
the confining potential (integrable or chaotic). We find that the peak bunch- 
ing effect seen in the experiments is reproduced in the theoretical models 
under certain special conditions. Although the independent particle model 
fails, in general, to explain quantitatively the short-range part of the peak 
height correlations observed experimentally, we argue that it allows for an 
understanding of the long-range part. 
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I. INTRODUCTION 



Since the first electronic transport measurements in semiconductor quantum dots in the 
Coulomb blockade regime were made,! a significant amount of data has been accumulated. 
The appropriate theoretical model used in understanding this data depends on the size of 
the electronic island. Two situations arise. When dealing with small or very small systems, 
where the number of electrons N in the Coulomb island is of the order of 10 (or a maximum 
of 10 2 ), the self-consistent approach is successfully In this regime the environment generates 
a smooth, nearly parabolic potential and most of the quantum dot properties can be under- 
stood in terms of an unsophisticated self-consistent electronic structure calculation. On the 
other hand, the most adequate theoretical description for data taken from larger systems, 
with iV of order 10 3 (or at least larger than 10 2 ), is statistical. In distinction to the first case, 
for large dots the confining potential is no longer expected to be parabolic up to energies 
close to the Fermi surface. In fact, in large gated structures, the dot geometrical shape 
can be substantially modified by adjusting voltages. Within the single-particle picture, it is 
believed that the complexity of the confining potential yields a chaotic motion in the clas- 
sical limit, thus justifying a quantum mechanical modeling based on a statistical theory.&i 
Further arguments for such approach are provided by the influence of weak disorder and the 
fact (not always justifiable) that many-body correlations should bring only small corrections 
to the usual single-particle description of transport. 

The single-particle statistical theory! - ! is widely accepted because of its early success in 
predicting the distribution of peak heights in the single-channel tunneling regime, later ex- 
perimentally confirmed.Hi In the regime where the mean level spacing A is much larger than 
fcsT, the predicted distribution of peak heights is strongly non-Gaussian and clearly peaked 
towards zero values. While the single-particle plus constant interaction theory captures most 
qualitative aspects observed in the experiments, some quantitative aspects defy comprehen- 
sion. Perhaps the most celebrated one is the absence of Wigner-Dyson fluctuations in the 
conductance peak spacingsJ^nL3 

There is also an important and rather prominent feature of the experimental data which 
is less discussed in the literature and remains unexplained: the notable clustering of peaks 
with similar heights. That is, in most experiments (an important exception is Ref. §), by 
sweeping the gate voltage, large and small peaks tend to appear in bunches, indicating the 
existence of peak to peak correlation. This is at odds with the statistical single-particle 
theory at low temperatures, which associates to each peak an independent single-particle 
eigenstate. A recent experimental showed yet another puzzling feature: the correlation 
tends to increase with temperature and the slope depends on the dot area (for a small dot 
it saturated at UbT w A/2.) 

Any successful theoretical model describing all these issues has to take into account that 
each conductance peak corresponds to a system containing a different number of particles. 
As electrons are added to the quantum dot the gate voltage configuration varies and the 
confining potential changes accordingly.0 Furthermore, and probably even more important, 
the additional particles change the electron mean field. Although not the original idea 
addressed in Ref. [14], we can parameterize the total change in the single-particle energies 
and wave functions by moving from one peak to the next through a parametric variation 
5X which takes into account changes of the external and mean field potentials. We can call 
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it the mean field "deformation". Although difficult to quantify, 5X is obviously nonzero. 

In this paper we concentrate on the statistical measures of the conductance peak heights, 
namely their distribution and autocorrelation function. We use the constant interaction plus 
independent particle model to discuss the role of "deformations" and temperature. We show 
that for ksT < A the theoretical peak height distribution agrees with the experimental data 
for systems with chaotic underlying dynamics, irrespective of the deformation parameter. 
We also show a very reasonable agreement between chaotic and some integrable models 
for ksT > 0.25A. As for the peak height autocorrelation function, we could not find any 
robust and generic mechanism to explain the strong peak to peak correlations seen in some 
experiments. We were only able to find such correlations for special choices of lead positions. 

The structure of this paper is as follows. In section |T| we begin with a short review of the 
constant interaction model, which is used throughout this work, centered on a discussion 
about finite temperature corrections and parametric variations. In section |TXT] we study the 
combined effect of temperature and deformations in a random matrix model. In section 



IV we investigate dynamical mechanisms to explain the observed correlations between peak 



heights. The proposed mechanisms are examined by studying a suitable dynamical model. 
The influence of data set sizes on the evaluation of peak height correlations is critically 
discussed in section [V]. Having analyzed only chaotic dynamical systems, we investigate 
in section |V| whether the experiments carried out so far indeed observed fingerprints of 
chaos. Insight on this question is obtained by studying the peak height distribution for a 



few integrable models. Finally, we summarize our conclusions in section VTJ 



II. THE SINGLE-PARTICLE PLUS CONSTANT INTERACTION MODEL FOR 

QUANTUM DOTS 

The theoretical description used in this study is based on the so-called independent 
particle plus "constant interaction" (CI) model,0 where the electron-electron interaction is 
taken into account through a fixed (capacitive) charging energy term in the Hamiltonian: 



N 



i=i 



u 



+ -N(N-l) + e V V g N, (1) 



where N is the total number of electrons with effective mass m* in the dot, U is the charging 
energy, and V g is some tunable gate voltage. The coefficient 77 is a function of the capacitance 
matrix elements of the dot .El The single particle potential in Eq. (|l|) accounts only for the 
electrostatic confinement and background disorder. 

When the dot is in equilibrium with the reservoirs represented by the leads, its total 
energy is minimal for a certain value of N. The precise value of N depends on the dot 
confining geometry, chemical potential (Fermi energy) in the reservoirs, and gate voltage V g . 
The electronic many-body ground state wave function is given by the Slater determinant of 
all bounded single-particle eigenstates with energies smaller than the Fermi energy. In its 
crudest version, the single-particle states in the CI model are determined by the bare electron 
confining potential alone. It has been arguedEHlil that many-body effects are important 
and one should instead get the single-particle states from a mean-field approximation, such 
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as Hartree-Fock. We will not attempt to answer this question in this work, but rather 
comprehensively study most features that may be traced to one-body effects. 

Both single-particle eigenenergies and eigenf unctions influence the electronic transport 
through the quantum dot. The wave functions, in particular, enter into the standard 
formulationi through the partial decay widths T r ' 1 which quantify the connectivity between 
the eigenstate in the dot and the electronic states in the leads (r and I stand for right and 
left leads, respectively). In the regime of thermally broadened resonances and large charging 
energy, r r,/ -C &;#T -C U, linear response theory yields the following expression for the the 
two-point conductance of a dot (weakly) coupled to reservoirs! 

2 oo "pZ "pr 

x [1 - f{SE -e N + e p )) 

xf(5E-TS(N)+TS(N - 1)), (2) 

with 5E = en + U(N — 1) — er]V g — Ef- The sum runs over all single-particle states p 
with energies e p and partial decay widths T r p and T l p . The latter describe the electron 
tunneling out of the dot and into the reservoirs, or vice-versa. Thus, the total tunneling 
rate is hY p = h(T l p + T p ). F eq (E p \N) is the canonical occupation probability of the level p 
for a dot containing N electrons in thermal equilibrium. / is the Fermi-Dirac distribution, 
f(s) = 1/[1 + exp(e//€BT)], S is the entropy of the dot, and Ef is the Fermi energy in the 
reservoirs. We will specialize our discussion to spin split (nondegenerate) energy levels. 

In experiments, sequences of conductance peaks are observed by varying the external 
gate voltage V g , thus changing the electrostatic energy in the dot. Within the CI model, 
in the single-level tunneling regime, where r r '' <C k^T <C A, and as the gate potential is 
swept, conductance peaks occur whenever the resonance condition 

E(N + 1) - E(N) = £ F (3) 

is satisfied.il Here E(N) is the ground state energy of the dot containing N electrons. Within 
the constant interaction model, Eq. (§) can be rewritten as 

e N + (N-l)U-er}V g = e F , (4) 

i.e., SE = 0. We remark that, although correct at the level of the Hartree-Fock approxima- 
tion (with AT ^> 1), it is still a matter of debate in the literature whether Eq. @ is generally 
valid. t3c3 The existence of many-body correlation effects may lead to a breakdown of Koop- 
mans' theorem and a different resonance condition. Furthermore, the strong assumption of 
a constant U prevents an adequate description of fluctuations in peak position. Neverthe- 
less, these restrictions shall not affect our discussion on the statistics of peak heights in the 
single-particle approximation. 

In order to make the notation more compact, let us rescale the conductance to a dimen- 
sionless form, 

r- e2 < r > ™ 
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Moreover, for the sake of simplicity, we consider that the average decay widths (taken 
over many resonances) are equal: (r r ) = (r ) = (r). In the temperature regime where 
(T) <C ksT <C A, the resonance condition, combined with Eq. (0), results in the well-known 
expression for the single-level conductance peak height, 



rir 



r 



a 



<'- g p (r) T l +T r ' 



Assuming that the electronic dynamics inside the dot is classically chaotic, Jalabert and 
collaborators^ used a random matrix model to calculate the distribution of partial decay 
widths in this temperature regime, in the presence (J3 = 2) or not {(3 = 1) of an applied 
magnetic field. They obtained 



PpO"®) = T^n ( TTT^Tm ) exp ( 1 , (7) 




where Cp is a normalization constant that can be expressed in terms of the gamma function: 
= (3/T(f3/2). By further assuming that the contacts are sufficiently far from each other so 
that r r and T l are uncorrelated, expressions were derived for the distribution of dimensionless 
peak heights, namely, 




(8) 
and 

P p=2 {a) = 4a[K (2a) + K^a)^ 2 ", (9) 

where Kq{x) and K\(x) are modified Bessel functions. These predictions were later experi- 
mentally confirmed by two independent experiments.!'! Equations (H) and (H) will serve us 
as benchmarks in later sections. 

As mentioned in the introduction, the early success of the random matrix hypothesis!'! 
in explaining statistical properties of conductance peaks in quantum dots seemed somewhat 
inconsistent. The experimental confirmation of theoretical prediction of the distribution 
P{q) contrasts with the observation of strong correlations between successive peak heights 
in most data sets.!lll Recall that standard random matrix theory (RMT)fH leads to uncor- 
related eigenfunctions, giving conductance peaks with a similar property. However, the fact 
that the condition ksT <C A is hardly achieved in the early experiments suggests that the 
discrepancy could be related to a finite temperature effect. The inclusion of finite temper- 
ature in RMT was carried out by Alhassid and collaborators^ by computing the canonical 
occupation factors F eq (e p \N) in Eq. @. Already for k B T w 0.3A they observed significant 
temperature corrections in P(a), as compared to Eq. ([]). The corrections, however, fail 
to quantitatively reproduce the available experimental datajll and the discrepancy grows 
larger as temperatures goes beyond A. We will return to this point later. 

The canonical average necessary to calculate F eq (e p \N) can be taken either by explicit 
summation of all possible single-particle configurations, as recently done in Ref. |25|, or in a 
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quadrature approximation for the canonical partition function, as in Refs. [24| , p6[ We adopt 
the latter approach in what follows. The procedure consists in expressing the canonical 
partition function as 

E tot (N)/k B T L AH-f , |A 

Z N = - £ I] + (10) 

where 0, = 27rj/L, £ tot (iV) = and 

1, (i>n 

with L standing for the number of states taken into account (we are interested in the 
case where N ^> L ^> 1). This form of organizing the sum, involving only exponentially 
decaying terms, together with the choice A = (ejy + £jv+i)/2, guarantees fast convergence. 
The free energy and the entropy are straightforwardly evaluated from Eq. (|10]) recalling 
that F(N) = -k B TlnZ N = E tot (N) - TS{N). With the quadrature approximation, the 
calculation of the canonical occupation number is carried according to 

L 

2 / ■ , lA-ei 



p -S(N)/k B L N +\ 

F eq (e P \N) = £ n (1 

j=ll =N -L 



x (l + e-^+^i^ . (12) 

Since we are only interested in the statistics of conductance peak maxima and not of 
spacings, we introduce in Eq. @ the resonance condition for k B T <^ A [Eq. @)]. In the 
limit of strong Coulomb blockade ksT e 2 /C, this yieldsll 

4 N +it Y l Y r 
x[l - f(e p -e N )) 

xf(-TS(N)+TS(N -1)). (13) 



One further remark about the implementation of Eq. (|10| ) is in order. In our numerical 
simulations we varied the number of eigenstates L summed over depending on the temper- 
ature range. For k B T < 0.1 A, we found that as few as 4 were sufficient to achieve good 
convergence; for higher temperatures, k B T > A, we used as many as 20. 

In the experimental setups where plunger gates are used to define both shape and depth 
of the confiningnotential, one expects that variations in the potential V g will continuously 
deform the dotE3 This situation is met in lateral quantum dots, where most of the data 
takenlilHl (an important exception is Ref. |28| ). Furthermore, as pointed out in 



were 



the introduction, as one adds electrons to the dot the mean field also changes. Thus, 
to establish theoretically the direct relation between V g and the dot electronic mean field 
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it is necessary to know the experimental set up in detail and perform rather sophisticated 
numerical simulations. While such study is desirable, it is not crucial for our analysis and will 
not be pursued here. We therefore simplify the problem parameterizing the "deformations" 
by a generic variable X on which eigenvalues and eigenfunctions will depend. In summary, 
the key point behind this hypothesis is that by adding one extra electron to the dot the mean 
electrostatic potential changes parametrically, on average, by 5XM In this way, the iVth 
conductance peak involves a sum over eigenergies and partial widths {e p (X), T r p (X), T l p (X)}, 
while the next peak in the sequence will be depend on {e p (X+SX), T p (X+SX), T l p (X+5X)}. 
It is convenient to measure the deformations in units of the typical distance between energy 
level anti-crossings, namely, the inverse root mean square derivative of the energy levels with 
respect to deformations: cB 

-1/2 

(14) 




X C = L 

(Here (■ • •) stands for an average over states p and over the parameter X). When the 
addition of a single extra electron is enough to scramble strongly energy levels and wave 
functions, Sx = 5X/X C > 1. In this circumstances, we expect each new conductance peak to 
signal a different and independent realization of the effective confining potential in the dot. 
Thus, in the regime of strong deformations, peaks heights should be fairly uncorrelated, in 
contrast with what happens as temperature is increased. A systematic study of temperature 
and shape deformation combined is carried out below. 



III. THE EFFECT OF PARAMETRIC VARIATIONS ON THERMAL AVERAGES 

In this section we investigate the effect of parametric Hamiltonian variations on thermal 
averages. More specifically, we are interested in answering two questions: (a) Are the 
correlations between consecutive peak heights sensitive to shape deformations? (b) What 
is the effect of parametric variations on the distribution of peak heights? We assume that 
the electronic underlying classical dynamics is fully chaotic and, therefore, its quantum 
Hamiltonian can be modeled by a random matrix. In doing so, we restrict our analysis to 
universal aspects; we leave the study of some specific dynamical models and nonuniversal 
features to the following sections. Since a fully analytical treatment of the problem is not 
available, we implement the numerical procedure discussed in the Sec. || to evaluate the 
canonical averages from a given sequence of energy levels and tunneling rates. 

The set of eigenenergies and eigenfunctions of a random parametric Hamiltonian can be 
conveniently obtained from a Hermitian matrix of the form 

H(X) = H lC osX + H 2 smX, (15) 

where Hi and H 2 are independent N x N matrices (N = 500 in our calculations) belonging 
to the proper Gaussian ensemble. In our numerical calculations, the orthogonal ensemble 
(GOE, (3 = 1) was used to model the situation when no magnetic field is present (time- 
reversal symmetric systems). For the case with non-zero magnetic field and broken time- 
reversal symmetry, we used instead the unitary ensemble (GUE, (3 = 2). For each realization 
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of H\ and H 2 and a given value of the parameter X, the matrix H(X) was diagonalized 
through standard methods. Only the 100 central eigenvalues were kept after each diago- 
nalization. This was done in order to avoid large density variations in the spectra and the 
consequent need for unfolding, as well as the inclusion of localized eigenstates from the tails 
of the band. The resulting eigenfunctions {ip p (k;X)} were used to generate the right and 
left tunneling rates according to the point contact approximation T r p (X) oc \ip p (k = 1;X)\ 2 
and Tp(X) cx \ip p (k = N; X)\ 2 . The energy eigenvalues {e p (X)} were used in the quadrature 
approximation for the canonical partition function given by Eq. fllPl). Thus, for a given real- 
ization of Hi and H 2 , a certain deformation 5X, and temperature T, we simulated sequences 
of dimensionless conductance peak heights as given by Eq. (|T3|). 

It is noteworthy that the calculation shown in Ref. [24], although similar in formulation 
to ours, relied only on the randomness of the tunneling rates and neglected fluctuations 
in the eigenenergies by adopting a picket-fence approximation for the spectrum. Here, 
instead, we use eigenenergies fully consistent with the Wigner-Dyson statistics predicted 
from RMT. For the unitary case, due to the strong level repulsion the picket-fence is a 
reasonable first approximation for the spectrum and we expect our results to be similar to 
those of Ref. Indeed we observe only small, 10% differences in peak height correlation 
functions (see below). For the orthogonal case, however, level repulsion is weaker and to 
obtain accurate results it is recommended to use the exact spectrum. This issue becomes 



very important when calculating F eq (e p \N) for dynamical systems. In Sec. [IV| we shall 
resume this discussion. 

The results presented here involved data obtained from 50 independent realizations of 
the matrices Hi and H 2 . For each of these realizations, 150 values of X were used within 
the interval [0, 7r/2] to generate peak sequences related to different "shape" deformations. 
Additional statistics was obtained by varying the "chemical potential" (i.e., the initial num- 
ber of electrons in the dot) and the "original shape" (i.e, the initial X point), taking care 
to avoid introducing spurious correlations between distinct peak sequences. 

As an illustration, three typical sequences of conductance peak maxima are shown in 
Fig. [l] for the orthogonal ensemble with some arbitrary deformation. Notice that, as tem- 
perature increases, the overall modulation becomes more pronounced, marking the existence 
of large peak-to-peak correlations. The clustering of peaks with similar heights is character- 
istic of the dominance of a single well-connected eigenstate over a wide energy range. 

The most characteristic feature of the thermal average is that the occurrence of small 
peaks becomes unlikely, even for temperatures which are smaller than A. This can be 



observed in the curves presented in Ref. |24| for the distribution of peak heights, although the 
effect is not always quantitatively confirmed by the experiments. For instance, in the case 
of preserved time-reversal symmetry the absence of small peaks should be very pronounced, 
in clear contradiction with the currently available experimental data.ill 

To start a quantitative analysis, let us first concentrate on the correlation of peak heights 
between neighboring conductance peaks. For this purpose, we have calculated the correlation 
function 

_ (5a N+n 5a N ) ^ 
where 5a n = «jv-(«jv) represents the deviation from average height of the iVth peak. The 
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averages (• • •) are taken over different realizations of Hi and H2, as well as over N. The 
results for each Gaussian ensemble at different temperatures and deformations are shown in 
Figs. |2] and [| Notice that even small parametric changes in the Hamiltonian can rapidly 
destroy the correlation of peak heights obtained from thermal averaging. This effect can 
be made more quantitative by calculating the half width at half maximum (HWHM) of the 
correlator c(n) as a function of temperature (see Fig. f|). 

Small deformations 5x exclusively due to shape variations are reported in a recent 
experiment!^ where, in addition, a theoretical curve solely based on "static" Gaussian en- 
sembles (5x = 0) is used to explain the data. Although the parametric deformation for this 
experimental set up is relatively small, one should not disregard its importance: For 5x ~ 0.3 
we already observe in Fig. § a much slower increase in correlation versus temperature than 



for 5x ~ 0, in agreement with the data of Ref. [13]. A fully consistent, parametric RMT 



for the peak height fluctuations yields a smaller correlation length than what is predicted 
by the "static" Gaussian ensembles. This is a clear indication that the theoretical curve 
shown together with the correlation length data in Ref. [13] should be taken with reservation. 
Nevertheless, the inclusion of parametric deformations alone does not solve the discrepancy 
between the RMT predictions and the experimental data at small temperatures. 

Another important statistical measure of the peak height fluctuation is its distribution. In 
Figs. |5| and |6| we show the distributions for both GOE and GUE cases. The inset illustrates 
the dependence on deformations for a given temperature. We observe that, contrary to 
the correlation lengths of peak heights, the distributions are not affected by parametric 
deformations, but do depend strongly on temperature. For ksT 3> A the distribution moves 
towards a Gaussian shape. It is important to notice that RMT captures the qualitative 
aspects of the temperature dependence observed experimentally,!!! but fails quantitatively: 
The data shows a smaller fluctuation in peak heights at high temperatures (ksT > A) than 
the theory predicts. This could be interpreted as an indication that for such temperatures 
many-body (or dephasing) effects become important. El 

We conclude that a theory based solely on random matrices with a proper thermal average 
is not capable of explaining simultaneously the peak height distribution and correlations. 
Even before attempting to improve the theory by including many-body interaction effects, 
we shall investigate whether nonuniversal aspects of the underlying single-particle dynamics 
also play an important role and lead to additional peak height correlation mechanisms. 



IV. DYNAMICAL CORRELATIONS IN CHAOTIC SYSTEMS 

In this section we investigate how single-particle dynamical correlations statistically in- 
fluence the conductance peak heights. We first examine how universal chaotic wave function 
correlations!^!! manifest on the conductance autocorrelation function. These correlations 
go beyond RMT theory, but are universal in the sense that they arise for any chaotic quan- 
tum state with finite wavelength. We show analytically that at T = the effect is small, 
but is enhanced at higher temperatures. We also investigate if peak height bunching can 
be explained by eigenstates localized in coordinate space. Such localization corresponds to 
wave functions with amplitudes concentrated along certain periodic orbits of the underlying 
classical Hamiltonian. This non- universal, spatial structure is known in the quantum chaos 
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literature as a "scar"a of the classical dynamics on the wave functions. 

We illustrate our investigation by studying numerically a suitable dynamical model. 
Namely, we calculate P(a) and c(n) for the smooth stadium potential as a function of 
T and deformation. We also model the coupling between the dot and the leads in two 
distinct ways. Most of the results presented below assume point contacts, in which case 
Tp oc |^ p (r c )| 2 , where ip p (r c ) is the p-th dot eigenfunction, evaluated at some point r c in 
the interface region between dot and lead. This description implies that, close to the dot, 
the leads should be narrow and carry a single propagating channel. Alternatively, we use 
the standard .R-matrix formalism to calculate the decay width for extended leads by taking 
the overlap integral over the contact area: T c p oc | / dsx c {r)tp p (r c )\ 2 , where \ c is the channel 
c wave function. Both models assume that the barrier penetration factors are small, such 
that (r) <C A. Furthermore, it is assumed that the penetration factors vary smoothly with 
energy and can thus be incorporated in (r) without affecting fluctuations. 



A. Corrections to the conductance fluctuations due to wave function correlations 

A recent work by Narimanov and collaborators!! explored the idea that the short time 
dynamics can influence peak height conductance fluctuations. Using the semiclassical theory 
they derived a correction to the conductance in the absence of magnetic field and ksT = 
which provides a particularly strong effect on the tails of the autocorrelation function. One 
of the important constraints in Ref. [35] is that the results are obtained for symmetrically 
placed leads, such that T r = T . Our study addresses a more generic situation and has a 
different starting point than theirs. 

In 1977, BerryS suggested that wave functions of billiards with a chaotic underlying 
dynamics can be statistically described by a random superposition of plane waves with a 
fixed energy. As a result, he found that chaotic wave functions at energy E display universal 
spatial correlations. In two dimensions the result is 

C(n, r 2 ; E) = A- l Uk\n - r 2 \), (17) 

where A is the billiard area, hk = (2m*£) 1/2 , and J (x) is the ordinary Bessel function. 
Numerical verifications of Eq. ( fl7|) usually use an arbitrary wave function, fix one coordinate 
and average over all directions to obtain C(5r,E), with 5r = \r± — r 2 \. Alternatively, and 
closely related to this study, C(5r, E) can be obtained by keeping r\ and r 2 fixed and taking 
the average over wave functions corresponding to eigenenergies close to E (see for instance, 
Ref. ^ and references therein). 

Since wave functions are spatially correlated, so are T l and T r . The joint probability 
distribution of r r and T l has already been obtained in the the orthogonal case (/? = 1) by 
the supersymmetric technique. El SrednickH later derived the same quantity by elegantly 
extending the formalism presented in Ref. We can summarize the result in the following 
way: Within the point contact model, defining v r ^ = T r ^'/{T), one findsffl 

Pp=i(v r ,v l ) = — J2N1/2 ,.. r .. m/2 exp 



2tt(1 - p)V 2 (v r v i y/ 2 



v r + v l 
2(1 -f 2 



V v 



x cosh \ [_p \ » ( 18 ) 
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where / = AC(6r, E) = Jo(k5r) in two dimensions. For / = one immediately realizes 
that the standard product of Porter-Thomas distributions is recovered. Unfortunately, we 
did not succeed in computing analytically the conductance peak height distribution starting 
from Eq. (f[8D. However, it is rather straightforward to express the average conductance as 
a function of / in terms of special functions. To leading order in powers of f 2 one finds 

(a(k5r)}, =1 = \{l + \[Uk5r)f} + 0(f). (19) 

For kbr > 1 we recover the average predicted by Eq. (|8p, as we should. In a typical 
experiment, by adding electrons and changing kdr, (a(k5r)) becomes a slowly oscillatory 
function on the energy scale of peak spacings. These long range oscillations in (a) do not 
appreciably change the autocorrelation function c(n) for small n, being only pronounced at 
the tails. These statement is in agreement with the results shown in Ref. By identifying 
k = kp (the Fermi wave number) and taking a typical experimental value for kpSr, the 
amplitude of the modulation will be of the order 1% of the standard RMT value for (a). 
This means that the maximum anticorrelation in c(n), at n = n anti , is also about 1% of c(0). 

A similar calculation can be made for the case of unitary symmetry ((3 = 2). Also in 
this case the joint probability distribution has already been derivedE3ta and reads 

pp=2(v ,v) = 7371 ex p [~T~r ) ( 1-/2 J > ( 2 °) 

where Iq(x) is the modified Bessel function. Here, again, it is difficult to obtain the conduc- 
tance peak height distribution P(ct) in a closed form. The averaged (a) as a function of /, 
on the other hand, can be found without much effort. The result in ascending powers of f 2 
reads 

(a(k5r)) p=2 = i {l + ^[Jo(Mr)] 2 } + 0(f), (21) 
which is an even smaller correction than for the orthogonal case. 



In distinction to the findings of Ref. |35], such spatial correlation effects will be hardly 
observable at zero temperature. We attribute the difference between our results and theirs 
mainly to their stringent constraint of having an exact symmetry in the wave functions near 
the tunneling region. Nonetheless, we anticipate that the long-range autocorrelation oscil- 
lations become more pronounced with increasing temperature, as we show in our numerical 
calculations. For future comparison with experiments, it is important to notice that the 
presence of an obstacle between r% and t*2 just requires the replacement of Sr by the length 
of the shortest classical path connecting 7*1 to r 2 HI Before concluding, an additional com- 
ment is in order. Until the present, Coulomb blockade experiments have offered very few 
fingerprints of the underlying classical electronic dynamics. An experimental check of the ef- 
fect proposed in Ref. exploring the symmetries of the wave function, would offer a clear 
indication of the validity range of the single-particle approximation and the semiclassical 
analysis. 
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B. Non-universal dynamical correlations due to scars 



Strong wave function correlations due to classical structures in phase (and coordi- 
nate) space, such as scars,@ can also lead to a non-universal behavior in the peak height 
fluctuations.il The analysis presented in this section is motivated by a recent work of Hack- 
enbroich and collaborators!! which explored the idea of shape deformation to explain correla- 
tions between heights of neighboring peaks. These authors proposed a correlation mechanism 
tailored for integrable systems. In what follows, we recast and reinterpret their mechanism 
for chaotic systems with strong scarring. Our results for this case are qualitatively distinct 
from those of Ref. |27|. 

For an integrable system it is well known that the structure of a given wave function 
does not change upon a parametric variation of the Hamiltonian. The reason is quite simple. 
The wave functions of integrable systems are strongly concentrated around classical tori, and 
the latter are only smoothly distorted by a parametric variation of the Hamiltonian. This 
behavior is generic, provided the system stays integrable upon the parametric variation. 
However, as the parameter X is varied, levels cross each other. Since they correspond to 
different sets of quantum numbers, the typical crossing distance is much smaller than the 
distance necessary to appreciably change the wave functions. Let us assume that in certain 
system there are levels very weakly influenced by changes in X. We call them "horizontal" 
levels h. (Any h level certainly has a quite marked wave function structure with respect 
to the system geometry.) "Horizontal" levels may cross generic levels, as the parameter X 
is varied. Let us assume now that the conductance peak of a dot containing N electrons 
is dominated by a certain h level for some value of X. If the typical 5X caused by the 
addition of one extra electron to the dot is of the order of the average parametric distance 
between crossings, the h level will always stay very close to the Fermi energy. In this way, 
the sequence of conductance peaks are dominated by the very same "horizontal" level and 
are thus expected to show correlations in heights. This mechanism!! is robust as long as 
the system dynamics is close to integrable. The restrictive condition of having exactly one 
anticrossing for 5X can be relaxed if temperature is included.il 

In the case of chaotic systems, wave functions typically decorrelate after one level 
crossing.!! Deviations of such universal behavior indicate the presence of scarred eigen- 
states. For such wave functions the mechanism described in the previous paragraph is 
applicable. However, one should be aware that if peak height correlations are explained by 
a non-universal feature they are not ubiquitous. 

In what follows we shall exemplify such mechanism and investigate how often these 
situation occurs. Moreover, by changing shape, position of the contacts, and temperature, 
we nearly exhaust all sources of peak height correlations within the single-particle scenario. 
The results are compared with those of RMT and the available experimental data. 



C. Dynamical model: the smooth stadium potential 

In order to investigate both aforementioned mechanisms, we studied the two-dimensional 
model Hamiltonian H = p 2 /2 + V(r) (we take the electron mass m* = 1), with the potential 
given byE3 
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(22) 



The exponent n sets the steepness of the confining potential. The model is flexible enough 
to allow for different classical behaviors. For instance, if we take n — 1 the system is 
integrable. As we increase the value of n our model Hamiltonian becomes very similar to 
the well-known stadium billiard, one of the paradigms of classical chaotic systems. For a = 1 
the equipotential V(x, y) — 1 corresponds to the border of the hard walled stadium billiard 
with unit radius and unit rectangular length. 

We model changes in the single-particle potential by varying a, the length of the rect- 
angular part of the well. The classical motion is chaotic if a is varied in the interval range 
1.0 < a < 1.25, provided we choose n = 2 and stay in the energy window around E — 1. 
Indeed, for this energy and parameter values the classical phase space is mostly chaotic, as 
can be seen in Fig. [7], though small remnants of integrability still exist. 

Finally, to assure that the single-particle levels we analyze correspond to a classically 
chaotic motion, we fix ft in order to have a window of energy around E < 1 containing 
typically 100 eigenstates of our Hamiltonian. The quantum eigenstates are obtained by 
numerical diagonalization. The problem is solved by using as basis states the eigenfunctions 
of a sufficiently large rectangular box. When the index n is an integer number, the calculation 
of the matrix elements is straightforward. The secular matrices are taken large enough such 
as to guarantee the convergence of all eigenvalues we analyze within 1% of the mean level 
spacing. This gives us confidence that the computed wave functions have converged. After 
ordering the eigenvalues in ascending energies, we consider only the eigenstates between the 
200th and the 300th levels. A representative region of the spectrum of the smooth stadium 
with n = 2, in the absence of an external magnetic field, is shown in Fig. §. 

Quantum calculations including a magnetic field B perpendicular to the stadium are also 
of simple implementation in this model. By the usual minimum ansatz p — > p + eA(r)/c 
with B = V x A we obtain a new Hamiltonian that can be diagonalized following the 
same numerical procedure. We found useful to work in the gauge where A = B/2[— c x (y — 
y ), c y (x— x )} with c x +c y = 2. Of course, the converged results do not depend on a particular 
choice of the gauge. However, the numerics can be made more efficient by choosing (xo,yo) 
close to the geometrical center of the stadium and c x L x = c y L y , where L x and L y are the 
lengths of the classically allowed region of the billiard at E = 1. With this choice, we 
obtain the largest ratio of number of converged eigenvalues to size of the secular matrix. We 
followed Ref. ET] to estimate the magnitude of the time-reversal breaking magnetic flux 



c 



threading the stadium. At the Fermi energy we obtain c /</>o ~ 2, where </>o is the unit flux 
quantum. When comparing the results of the dynamical model against the RMT unitary 
case, we fix the magnetic field to be BA = <f) c , where A is the classical allowed area for an 
electron with energy Bp. 

Although the system is chaotic, there are many very narrow avoided crossings, which are 
characterized by a gap e< A. These crossings occur for only a few eigenstates, indicating 
that the latter are very weakly coupled to others eigenstates. This non-universal feature is 
well-known in quantum chaos and is a fingerprint of scarred wave functions. Scarred eigen- 
states have the property of concentrating wave function amplitudes along certain classical 
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periodic orbits. Upon a parametric variation and after a level crossing, scarred eigenstates 
preserve their strong identity. In view of Fig. || we can apply the mechanism proposed for 
integrable systems in Ref. to explain strong peak-to-peak correlations. We identify the 
h levels with scarred "horizontal" states. The states that contribute to the conductance in 
the case where SX = X c and k B T <C A are indicated in Fig. |] by dots. The two sequences 
shown correspond to different Fermi energies. The wave functions related to the lower en- 
ergy sequence are concentrated along the family of periodic orbits that bounce from y = 
to y — 1 with v x = (the so called bouncing ball modes). The corresponding conductance 
peak heights show very large peak-to-peak correlations when the leads are placed such as 
to intersect these bouncing trajectories. The result is shown in Fig. ^. Notice that for 
small temperatures the mechanism holds even if the scarred state does not correspond to 
the Fermi level, but is its neighbor. This happens because scarred levels show very weak 
level repulsion and very often contribute significantly to the canonical occupation factor at 
k B T <C A. (Differences in c{n) due to spectral fluctuations were already discussed in pip. 

The proposed mechanism holds for chaotic systems if three main conditions are met. 
First, the tunneling state at the Fermi energy has to be localized in coordinate space. 
Second, SX has to be close to X c to keep this state close to ep. Third, the leads have to 
be in the region where the wave function has an enhanced amplitude. In practice, these 
conditions are obviously difficult to be satisfied simultaneously. However, one could imagine 
that a given sequence can have few correlated peaks dominated by a particular scarred state. 
When this state drifts away from the Fermi level it will eventually be replaced by another 
localized state. In this way a peak height sequence would have several regions where peaks 
are strongly correlated. This is the kind of effect that we actually tried to identify. Although 
we observed sequences like that in the numerical simulations, the proposed mechanism is not 
generically robust. When one averages over many sequences by varying the positions of the 
contacts inside the billiard, the Fermi level, as well as 5X C . The autocorrelation functions 
c(n) dependence on k B T/A does not differ much from the RMT result. In other words, 
sequences like the top one in Fig. § are statistically dominant. 

Summarizing: chaotic systems like the stadium can, in some cases, display strong peak- 
to-peak correlations. Such cases, however, are statistically rare and in our opinion this 
mechanism is not robust enough to explain the experimental data.Uil 

Let us now turn to the study of long range correlations. For this purpose we evaluated 
c(n) for the smooth stadium without magnetic field. Deformations which cause energy levels 
to change by A are very small at the classical level and do not influence our results. We 
therefore only consider the case of 5X = 0. The data set consists of 60 sequences of 100 
peaks each. The peaks in each sequence are normalized to unit average height. This is 
done to avoid the systematic errors due to the decrease of the mean conductance from one 
sequence to another (since each sequence corresponds to a different length a). A correlator 
is calculated for each sequence and an average over all sequences is taken. For k B T = no 



undershooting is observed in c(n), as predicted in Sec. [IV A[ The results for k B T = 0.5A are 



shown in Fig. ITUl We do not know how to estimate the magnitude of the oscillations in c(n) 



after the temperature average, since / enters in the expressions in a convoluted manner. In 
Fig. [n] we observe much larger oscillations than predicted by the theory at k B T = 0. We 
can check if the enhanced oscillations are related to the mechanism proposed in [IV A| , by 
verifying if n anti (the position of the anticorrelation maximum) is consistent with our theory. 
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Our data were obtained by sweeping kp and keeping the lead positions fixed. For spinless 
electrons in a billiard, k F « (4nN/A) 1/2 . To add one electron to the dot, we have to move 
up in energy by A in our model. This corresponds to change kp by lisjAkp. Using the 
parameters of our calculations, we then estimate the position of the maximal anticorrelation 
to be ra ant i ~ 8, which has to be contrasted with Fig. |1(]. A similar exercise was done in 
Ref. [35| using realistic values to estimate n ant i from the experimental data, obtaining values 
in reasonable agreement with Ref. O. In our work, anticorrelations in c(n) cannot be ruled 
out for small values of fcgT/A; nevertheless, they have to increase at higher temperatures 
and this is observed in our numerical simulations. We postpone the presentation of this 
result to the next section, where we also discuss the possible statistical problems inherent 
to such study. 



V. FINITE DATA SET AND THE PEAK HEIGHT AUTOCORRELATION 

FUNCTION 

The proper statistical treatment of experimental data is of great importance in any 
fluctuation analysis. In practical situations the number of Coulomb blockade peaks in a 
sequence goes from 50 to 100. Therefore, the problem of biases and errors due to finite sample 
size is one of the major difficulties in calculating peak height autocorrelation functions. If a 
reliable estimate of the effective number of independent points N e s can be made, the problem 
of making corrections and estimating errors is partially solved. 

Considerable effort was devoted to this type of problem in the context of compound 
nucleus reactions (see, for instance, Ref. ^ and references therein). In this case the observ- 
able of interest is the nuclear cross section, whose autocorrelation function is known to be 
Lorentzian. Let us translate this earlier estimate for the effective number of independent 
peaks into the language of Coulomb blockade. For this purpose, we assume momentarily that 
the peak height autocorrelation function has a Lorentzian shape. Doing so and introducing 
r = n/n c , we can write the variance of the average conductance peak height as 



var (a = var a 



2 tan l r ln(l + r 2 ) 

(23) 



where var(a) is the variance of the conductance peak height (either the "true" variance or, 
for an estimate, the ensemble variance), whereas a is the average over a sequence of n peaks 
with a correlation length n c . It is known that if an average of iV e g- independent data points 
is computed, the variance of this average is given by 

. . var(a) 

var(a) = — M. (24) 

By comparing the two last equations, we can express iV e fj in terms of n/n c . For values of 
n/n c ^> 1, A/gff is very well approximated by 

In , . 

iV eff « , 25 

■nrir. 
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which tells us that an effective independent point is found after an interval of irn c peaks. 
The statistical fluctuations are, as usual, given by l/^N e g. We expect our estimate of the 
number of effectively independent points to change little when a different functional form 
for the autocorrelation function is considered. The factor it will be replaced by another 
numerical factor of the same order of magnitude. 

We illustrate this discussion by considering peak height autocorrelation functions ob- 
tained from sequences of different lengths for the smooth stadium without magnetic field 
or deformation for ksT = 0.5A and 2A. As temperature increases, so does n c ; statistical 
fluctuations become large even for samples with many peaks. This is illustrated in Fig. [TTJ. 
At ksT = 2 A, n c = 5 and each sequence of 25 peaks has N e s = 1 and fluctuations prevail 
at n ^> 1. Similar results are obtained when a magnetic field is included, but are not shown 
here. 

Notice that the correlation length is more sensitive to small samples at high temperatures. 
For instance, at fc^T = 2A, n c decreases appreciably as we consider shorter sequences. 
However, at ksT = 0.5A the correlation length (defined as the HWHM) is stable. As 
a consequence, the curve n c (T) obtained from short sequences tends to saturate at high 



temperatures. This is similar to what is seen in Ref. [13| for the dot with A = 38u.eV (Dot 
2), the smallest one among the three dots tested. Our conclusion is that it is rather difficult 
to affirm unambiguously that the saturation of n c with temperature is a manifestation of 
decoherence. It may well be the effect of small statistics. 



VI. ARE CONDUCTANCE FLUCTUATIONS A FINGERPRINT OF CLASSICAL 

CHAOS? 

Any systematic search in the literature for single-particle models designed to explain 
the fluctuations of Coulomb blockade peak heights will reveal a strong bias towards the 
assumption of chaos. Very little attention has been devoted to integrable models, with, to 
the best of our knowledge, one exceptional This is quite understandable since it is believed 
that, due to small irregularities in the confining potential or the presence of weak impurities, 
the electronic dynamics tends to be chaotic. In fact, it is particularly difficult to set up an 
experimental realization of an integrable system. This statement is corroborated by the fact 
that there is only one kind of conductance experiment that shows a clear distinction between 
integrable and chaotic underlying dynamics.0 However, this study, carried out by Chang 
and collaborators, dealt with the weak localization peak in open quantum dots, i.e., not in 
the Coulomb blockade regime. 

Looking for clear fingerprints of integrable dynamics in the statistics of conductance 
in the Coulomb blockade regime is probably not a well posed task. Integrable systems 
are not universal, implying that each different confining potential will likely lead to very 
different conductance fluctuations. In this section we pursue the opposite strategy, asking 
if it is possible to tell from the conductance experimental data if the underlying dynamics 
is chaotic or not. To obtain some insight into this question we analyze two very simple 
systems, namely, the rectangular billiard (a paradigm of integrable motion) and the smooth 
stadium with n = 1. 

We model the coupling of the quantum dot to a given lead by a point contact at the 
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position r c , which implies that the decay width of the pth state is proportional to \ip p (r c )\ 2 . 
For billiards it is standard to solve the single-particle problem with Neumann boundary 
conditions and to put the point contacts at the boundaryaS For chaotic systems one obtains 
equivalent results choosing the point contacts inside the dot, thus probing the "bulk" of the 
single-particle wave functions. In this section we investigate both schemes. 

For the rectangular billiard we can only proceed analytically in the regime where T <C 
ksT <C A. The decay width of a state with quantum numbers (n,m) is given by T = 
(4/ L x L y ) sin 2 (k n x) sin 2 (k m y), up to a penetration factor. In the absence of an external 
magnetic field, in the "bulk" point contact model, P(T) can be calculated through the 
expression 



bulk 



\ rL x rL y . 4 9 

(r) = — — / dx dy 5{T - — — sin (k n x 
LtLti jo jo Ajrrijy 



xsin 2 (k m y)). (26) 



For point contacts, one immediately obtains (T) = (LxLy)' 1 . The integrals in Eq. fl26|) 
represent an average over the position of the point contact. One could object that this is 
not the kind of average taken in experiments. However, it turns out that by reducing the 
integration in Eq. (|26"D to a single period in x and y, the wave function quantum numbers 
drop out of the integral and P{T) becomes independent of n and m. In this sense, we claim 
that the average in Eq. (^) is equivalent to an average over a large number of wave functions 
with a fixed point contact position. 

Performing the integral in Eq. (|26| ) and writing the result in terms of (r), we arrive at 



i(Ulk (r) = — _ / 1 = (27) 



which can be written in a closed form in terms of a complete elliptic integral of the first kind 
(not shown here). 

In a similar manner, we can calculate the distribution of decay widths P(T) for the 
rectangular billiard with point contacts placed at its sides, obtaining 

p sidc (r) = - 1 (28) 
* Jr(2(T) - r) 



In Fig. [12] we compare the analytical results of P(T) for the rectangle with the Porter-Thomas 
distribution, which is the RMT prediction for the orthogonal symmetry. In addition, we also 
display the numerical results for the smooth integrable stadium potential (n = 1). 

The similarity between the bulk rectangle and the smooth stadium with the Porter- 
Thomas distribution (in the absence of magnetic field) can be understood by a simple argu- 
ment. Any wave function has a discrete number of points where its amplitude squared is a 
local maximum. It has also nodal lines where the amplitude is zero. By picking coordinates 
at random, the probability for obtaining a small amplitude is much bigger than for a large 
amplitude. It has to be emphasized, however, that neither the bulk rectangle nor the smooth 
stadium are in perfect agreement with the Porter-Thomas distribution. Although similar, 
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there are differences that can be better visualized in a log-log plot (not presented here), as 
shown before in the quantum chaos literature (see, for example, Ref. |43|) . 

The most interesting aspect of this analysis appears when we treat the distribution of 
conductance peak heights P(a). We find that P{ol) is relatively insensitive to P(T), as illus- 
trated in Fig. [H| Analytical calculations are rather difficult in this case, since, in general, 
we do not know how to take into account correlations between the decay widths r z and 
r r for integrable systems. Thus, the presented results were obtained numerically. Notice 
that although the histograms are somewhat different from the RMT result, they become 
indistinguishable when the usual experimental uncertainties (5-10%) are involved. (To con- 
struct the histogram shown in Fig. [13] we have used a very large sample of conductance peak 
heights.) The apparent similarity between the RMT prediction and the result for integrable 
systems becomes even stronger if we take k B T = 0.5A, as illustrated in Fig. Ob. 

From this analysis, we conclude that the universal aspect observed experimentally for 
P(a) is more robust than initially predicted. We believe that there is a large class of 
integrable systems which lead to conductance peak distributions similar to RMT. This result 
is not very intuitive, since for most observables integrable systems show larger fluctuations 
than chaotic ones.a We caution that our statement is not that any integrable system will 
display P(ot) similar to Eqs. (|§|) and (|9|). The literature already offers us a counter example: 
The conductance fluctuations observed in the model analyzed in Ref. 25 are certainly not 
compatible with RMT. 



VII. CONCLUSIONS 

In this work we aimed to study quantitatively how Coulomb blockade conductance peak 
statistics can be understood within the single-particle plus constant interaction model. We 
showed that the distribution of peak heights predicted by the RMT is very robust to de- 
formations in the dot effective potential. For k B T > 0.25A we found that there are exist 
integrable systems which give equivalent results to RMT. The curves derived numerically 
from RMT reproduce fairly the experimental results at low temperatures (for k B T < 0.5A), 
with increasing discrepancies as temperatures rises. 

We also showed that the peak height long-range correlations can be understood in terms 
universal wave function correlations. The short range correlations remain an open problem. 
Under certain circumstances, we were able to reproduce peak sequences with visible bunching 
of peak heights. Their correlation lengths depend on the specific portion of the spectrum 
considered and is thus nonuniversal. Unfortunately, such sequences are relatively rare, since 
they rely on the existence of well-connected scarred states. 

Starting from non-degenerate sequences of eigenenergies and tunneling rates, one can 
construct spin-degenerate data sets by simply duplicating the original sequences. In this 
case, at k B T/A — > the conductance peaks should appear in pairs of equal heights. More 
importantly, paired peaks would have a similar response to external perturbations, such as 
a magnetic field, even if their heights are not quite identical. Experimentally, however, this 
behavior is not observed (with the possible exception of Fig. 2 in Ref. [13]). In view of this 
fact, we have decided for restricting our analysis to nondegenerate single-particle levels. 

It is noteworthy that the experiment by Chang and collaborators! is the only one where 
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peak bunching is not visually strong. The main difference between the setup of Ref. |8| and 
others is the sample, which in this case has lower mobility. Disorder destroys the dynamical 
effects we propose. One could then speculate that the lack of peak height correlations in 
Ref. II and, perhaps, the departure from universality seen in other experiments depend on 
the strength of dynamical effects. From this perspective, our main conclusion is negative, 
since, we were not able to identify a robust mechanism to explain such phenomena based on 
this idea. 

The hypothesis that the ground state and the excited states are uncorrelated and follow 
RMT leads always to the suppression of small conductance peak heights upon thermal 
averaging. Even the inclusion of additional features, like deformations, does not change 
this conclusion. In view of the experimental evidences of an even larger suppression of the 
small peaks than predicted by RMT, we might conclude that the ground state and the 
first excited states have instead some correlation. The mechanism responsible for this could 
also be responsible for the correlations between ground states with consecutive number of 
electrons. Thus, peak bunching and the persistence of small peaks could be manifestations 
of the same mechanism. From our analysis, however, such mechanism seems to be beyond 
the mean-field scenario. In conclusion, we believe that the full quantitative understanding 
of these effects may lay outside the single-particle approximation. 
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FIG. 1. Sequences of dimensionless conductance peak maxima taken from the GOE simu- 
lations for a given set of eigenstates and no deformation at (a) fc#T = 0.1A, (b) ksT = A, and (c) 
k B T = 3A. 
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FIG. 2. Peak height correlation function for the GOE at temperatures ksT = 0.1A (circles), 
0.3A, 0.5A, A, 1.5A, 2A (stars) for different shape deformations: (a) 5x = 0, (b) 5x = 0.1, (b) 
Sx = 0.5, and (d) Sx = 1.0. 



23 




24 




25 




12 3 

a 



FIG. 5. Distribution of peak heights for the GOE without parametric distortion at different 
temperatures: ksT = 0.1A (solid line), 0.5A (dotted line), and 1.5A (dashed line). The inset 
shows the distribution at T = 0.5A obtained for three different deformations (8x = 0.1, 0.3, and 
1.0). The heights are rescaled to (a) = 1. 
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FIG. 6. The same as in Fig. | but for the GUE. 
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FIG. 7. Poincare surface of section for the smooth stadium billiard with n = 2, a = 1, and 
E = 0.75. The section is taken at x = a. 
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FIG. 8. Energy levels as a function of the length a for the smooth stadium with n = 2 and no 
magnetic field. The two nearly horizontal sequences of dots correspond to tunneling states at low 
temperatures regime and for Sx = 1. The lower energy sequence follows a scarred pattern, while 
the higher one does not. 
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FIG. 9. Sequence of peak heights for "bouncing ball" wave functions at fc^T = 0.25A. 
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FIG. 10. Peak height autocorrelation function c(n) for the smooth stadium (n = 2, no magnetic 
field) obtained from 60 sequences of 100 peaks each, going from ksT = 0.25A (circles) to 2A (stars) 
in increments of 0.25A. 
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FIG. 11. Peak height autocorrelation function c(n) for the smooth stadium obtained from 
different sets of sequences: 60 sequences of 100 peaks (crosses), 120 of 50 (squares), and 240 
sequences of 25 peaks (circles), (a) ksT = 0.5A; (b) ksT = 2A. 
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FIG. 12. -P(T) from RMT (solid line) compared with different integrable systems: bulk rectan- 
gle (dotted), side rectangle (dashed), and n = 1 smooth stadium (dash-dotted). The distributions 
are scaled to have (r) = 1. 
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FIG. 13. P{a) from RMT (solid line) compared with different integrable systems: bulk rect- 
angle (dotted), side rectangle (dashed), and n = 1 smooth stadium (dash-dotted) for (a) ksT = 
and (b) ksT = 1.0A. The distributions are scaled to have (a) = 1. 
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